Resonant addressing and manipulation of silicon vacancy qubits in silicon carbide 



D. Riedel\ F. Fuchs\ H. Kraus\ S. Vath\ A. Sperlich\ V. Dyakonov^'^, 
A. A. Soltamova^, P. G. Baranov^, V. A. Ilyin"^, and G. V. Astakhov^* 
^Experimental Physics VI, Julius- Maximilian University of Wiirzburg, 97074 Wiirzhurg, Germany 
^Bavarian Center for Applied Energy Research (ZAE Bayern), 97074 Wiirzhurg, Germany 
^loffe Physical- Technical Institute, 194021 St. Petersburg, Russia 
^Saint Petersburg Electrotechnical University, 197376 St. Petersburg, Russia 
(Dated: October 2, 2012) 

Several systems in the solid state have been suggested as promising candidates for spin-based 
quantum information processing. In spite of significant progress during the last decade, there is a 
search for new systems with higher potential [D. DiVincenzo, Nature Mat. 9, 468 (2010)]. We report 
that silicon vacancy defects in silicon carbide comprise the technological advantages of semiconductor 
quantum dots and the unique spin properties of the nitrogen-vacancy defects in diamond. Similar 
to atoms, the silicon vacancy qubits can be controlled under the double radio-optical resonance 
conditions, allowing for their selective addressing and manipulation. Furthermore, we reveal their 
long spin memory using pulsed magnetic resonance technique. All these results make silicon vacancy 
defects in silicon carbide very attractive for quantum applications. 



PACS numbers: 61.72.Hh, 71.55.-i, 76.70.Hb, 61.72.jd 

The double radio-optical resonance in atoms [1] con- 
stitutes the basis for a unprecedented level of coherent 
quantum control. Atomic time standards [2] and multi- 
qubit quantum logic gates [3] are among the most known 
examples. In the solid state, semiconductor quantum 
dots (QDs) and the nitrogen-vacancy (NV) defects in 
diamond, frequently referred to as artificial atoms, are 
considered as the most promising candidates for quan- 
tum information processing [4, 5]. Nevertheless, such a 
high degree of quantum control, as achieved in atoms, 
has not yet been demonstrated in these systems so far. 
Therefore, there is a search for quantum systems with 
even more potential [6]. 

Recently, intrinsic defects in silicon carbide (SiC) have 
been proposed as eligible candidates for qubits [7, 8] . In- 
deed, they reveal quantum spin coherence even at room 
temperature [9-11]. All of these experiments have been 
carried out under non-resonant optical excitation where 
all spins are controlled simultaneously. However, for spin- 
based information processing it is necessary to perform 
manipulations of selected spins, while the rest should re- 
main unaflFected. This demonstration in SiC is still an 
outstanding task. 

The selective spin control can be realized using a reso- 
nant optical excitation. As a rule, inhomogeneous broad- 
ening is much larger than the natural spectral linewidth, 
and such resonant addressing can be done on single cen- 
ters only. To avoid this problem, we applied a special 
procedure to "freeze" sificon vacancy (Vsi) defects dur- 
ing their growth, allowing to preserve a high homogeneity 
inherent to Lely crystals. This is confirmed by the ex- 
tremely sharp optical resonances in our samples. The 
spectral width of the Vsi absorption lines is several /ieV 
(ca. 1 GHz) , which comparable with that of a single QD 
or a single NV center in diamond. 

We then demonstrate the selective spin initialization 



and readout by tuning the laser wavelength together with 
the spin manipulation by means of electron spin reso- 
nance (ESR). Such a double radio-optical resonance con- 
trol indicates that the Vsi defects strongly interact with 
light and are well decoupled from lattice vibrations. The 
latter is also confirmed by the observation of a long spin- 
lattice relaxation time (Ti =0.1 ms). This can be po- 
tentially used for the physical implementation of scalable 
multi-qubit quantum logic gates in solid state devices 
[12] , vector magnetometry with nanometer resolution [13] 
and quantum telecommunications via spin-photon cou- 
pling [14]. 

SiC is a wide-gap semiconductor {Ef^~^'^^ = 3.05 eV) 
and possesses a spectrum of unique mechanical, electri- 
cal and thermophysical properties making it appropri- 
ate for many demanding applications [16]. The material 
properties of SiC, including intrinsic defects, are being 
investigated for decades. In particular, the Vsi defects 
have been identified in the 80's [17]. SiC exists in about 
250 crystalline forms, called polytypes, which are varia- 
tions of the same chemical compound that can be viewed 
as layers stacked in a certain sequence. Here we con- 
centrate on the polytype 6H-SiC [see Fig. 1(a)]. It has 
the stacking sequence ABCACB and is characterized by 
three nonequivalent crystallographic sites, one hexagonal 
(h) and two quasicubic (kl, k2). 

The investigated 6H-SiC samples of high crystalline 
quality have been grown by the modified Lely method. 
The high-temperature (2700° C) seedless crystal growth 
is driven by the temperature gradients within the cru- 
cible, resulting in a pressure gradient and thus, in a mass 
transport. The process is followed by the subsequent fast 
cooling, which "freezes" the defects within the lattice at 
low densities (for details see [15]). Taking into account 
the retrograde character of the nitrogen solubility with 
temperature, the doping level of the 6H-SiC crystals is 
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below 10^^ cm~^ and the compensation degree of nitro- 
gen donors is ca. 20%. The Vsi concentration hes below 
this value and we estimate that ca. 10^^ Vsi are probed 
in our experiments. 

In all experiments, the sample was mounted in a liq- 
uid helium flow cryostat with a temperature controller. 
PL was excited by a He-Ne laser and detected by a CCD 
camera linked to a 800-mm monochromator. For the res- 
onance excitation we used a tunable diode laser system 
with a linewidth below 1 MHz. The ESR experiments 
were performed in a home-modified X-band spectrome- 
ter, with the magnetic field applied parallel to the c axis 
of 6H-SiC. The magnetic field inside the spectrometer 
was modulated at a frequency of 100 kHz and the ESR 
detection was locked-in. 

Figure 1(b) presents a typical low-temperature (T = 
4K) photoluminescence (PL) spectrum of our 6H-SiC 
sample obtained with sub-bandgap excitation with a 
HeNe laser {E = 1.959 eV). The PL consists of sharp zero 
phonon lines (ZPLs) and their sideband phonon replicas. 
Three of these ZPLs - labeled as Vsi(h), Vsi(kl) and 
Vsi(k2) - originate from silicon vacancies at the corre- 
sponding crystallographic sites [18]. The origin of the 
UD-3 ZPL is still under debate [19], and we do not dis- 
cuss it in what follows. 

We now present light-enhanced electron spin resonance 
(ESR) experiments. It has been shown that the ground 
state of silicon vacancies is a high-spin state [20]. In 6H- 
SiC with hexagonal crystal structure the Vsi defects have 
the point group symmetry Cav, allowing for a zero-field 
spin splitting (A). The application of an external mag- 
netic field leads to a further splitting of the Vsi ground 
state. When the difference between spin-split sublevels 
meets the ESR frequency {vesr) the electron spin reso- 
nance occurs. The energy diagrams describing the optical 
transitions, the spin pumping scheme and the related RF 
transitions with respect to the multiplicity of the ground 
state are discussed in [21]. In case of the silicon vacancy, 
two ESR lines should appear for each Vsi site at mag- 
netic fields B- and 5+ [8, 21]. If the external magnetic 
field is applied parallel to the c axis of 6H-SiC, there is 
the following interconnection between i^esr and A 

hi^ESR = tA + gel^BB± . (1) 

Here, = 2.0 is the electron g-factor and 113 = 5.79 x 
10~^eV/T is the Bohr magneton. The amplitudes {A±) 
of the Vsi ESR lines depend on the population differ- 
ence between the particular spin sublevels involved. In 
the dark the difference is due to the Boltzmann factor, 
and the amplitudes A± are negligibly small. The optical 
excitation of Vsi defects and following relaxation pref- 
erentially pump the system into certain spin sublevels 
of the Vsi ground state. This results in light-enhanced 
ESR, as exactly observed in our experiments presented 
in Figs. 1(c) and (d). 
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FIG. 1: (color online), (a) The Vsi defects present at three 
nonequivalent crystallographic sites of the 6H-SiC lattice, 
one hexagonal (h) and two quasicubic (kl and k2). (b) 
PL spectrum obtained under excitation with a He-Ne laser 
[E = 1.959 eV). The ZPLs of the corresponding Vsi defects 
are labeled by arrows. The shaded area from 1.20 to 1.34 eV 
indicates the spectral range, where the ODMR of Fig. 4 is de- 
tected. (c,d) Light-enhanced electron spin resonance (ESR) 
under excitation into the Vsi(h) ZPL {Eh = 1.397 eV) and 
above-ZPL excitation {E — 1.412 eV), recorded at tempera- 
tures T = 50 and AK. The Vsi(h) ESR lines are observed at 
magnetic fields Bk- = 331.9 mT and Bh+ = 341.0 mT. The 
ESR experiments are performed in a X-band spectrometer 
(zy^5^ = 9.43451 GHz, B\\c). 



First, we discuss the high temperature data [Fig. 1(c)]. 
The three lines in the magnetic field range from 335.0 
to 337.5 mT are the well known ESR fingerprint of the 
nitrogen (N) donor in 6H-SiC [16], being independent of 
the optical excitation. The pair of outer lines at magnetic 
fields Bh_ = 331.9 mT and B^^ = 341.0 mT appears 
in the ESR spectrum under optical excitation into the 
Vsi(h) ZPL {Eh = 1.397eV). Using Eq. (1) we obtain 
the zero- field spin splitting in the Vsi(h) defect = 
0.527/ieV (127 MHz), which is in agreement with the 
earlier reported value [22]. 

At cryogenic temperatures, the Vsi(h) ESR lines are 
also observed under the above-ZPL excitation {E = 
1.412 eV), as shown in Fig. 1(d). In addition, the pair 
of Vsi(k2) resonances appears. Again using Eq. (1) we 
obtain the zero-field spin splitting in the Vsi(k2) defect 
Ak2 = 0.11 /ieV (27 MHz) [22]. 

Remarkably, under the optical excitation the B_ and 
5+-lines of Vsi (both h and k2) have opposite phase in 



3 



a 

(a)S^^ = 341 .0 mT 
H ; S ! ! — ( 


^ V3,(h) 


=331.9 mT 
r=4K 


— % 


b > 


@B^_^ 335.5 mT 



0.4 
^ 0.3 

I.. 

QQ 

I 0.1 
QQ 

P 0.0 

.0) 

i -0.1 
1-0.2 



0) -0.3 
DC 



-0.4 



c 




d 




-> 










Si " 






~ -> 

29o- 








- Si 

- 


■TP* - 




_ ^ 


m m 












-> 








-2x^^Si . 


- 




III: 





1.36 1.38 1.40 1.42 -400 -200 200 
Excitation energy, E (eV) Rel. exc. energy, E -E^ (laeV) 

FIG. 2: (color online), (a) The peak-to-peak values Ah± of 
the Vsi(h) ESR lines [see Fig. 1(c)] recorded in magnetic fields 
Bh- — 331.9 mT and Bh+ — 341.0 mT as a function of exci- 
tation energy, (b) The same as (a), but for Ak2- recorded in 
a magnetic field Bk2- = 335.5 mT. (c) A 2D color plot of the 
ESR signal as a function of the relative magnetic field B — Bu- 
and relative excitation energy E — Eu- Here, Bu- — 331.9 mT 
and Eh = 1.397 eV correspond to the ESR and optical reso- 
nance of the Vsi(h) defect, respectively, (d) A cross-section of 
c E — Eh = —7^ /ieV (as indicated by the vertical arrow). 



ESR spectra. This is a signature of the spin pumping 
under optical excitation, leading to the inverse popula- 
tion between spin sublevels. As a result, radiofrequency 
(RF) emission rather than absorption is detected for one 
of these ESR lines [21]. For Vsi(li) the RF emission is 
observed at the magnetic field 5/^-, as schematically in- 
dicated in Fig. 1(d). For convenience we assume positive 
value for the absorbing resonances > 0) and nega- 

tive value for the emitting resonances {A}^_ < 0). 

Figure 2(a) shows how Ah± depend on the excitation 
energy. For energies E < the spin pumping is ineffi- 
cient, for E > Eh it growths monotonically, presumably 
due to phonon-assisted processes and for E = E^ ^ very 
sharp optical resonance is detected. A similar behavior 
is observed for Ak2- with E = Ek2 [Fig. 2(b)]. 

We are now in the position to discuss the double radio- 
optical resonance of Vsi(h). The ESR signal as a function 
of the relative magnetic field B — Bh_ and relative ex- 
citation energy E — E^ is presented in Fig. 2(c). One 
can relate the ESR satellites at 5 — B^- = ±0.17 mT 
and ±0.34 mT to one and two nuclear spin-carrying iso- 
topes ^^Si (7 = 1/2) among the 12 next-nearest-neighbor 
silicon atoms [Fig. 2(d)]. The optical resonance reveals 
a complex structure consisting of a series of extremely 
sharp lines. It is more clearly seen in Fig. 3(a), where 
the Vsi(li) optical resonance from Fig. 2(a) is shown with 
higher resolution. Remarkably, the spectral width of an 
isolated line is about 2 /ieV [the inset of Fig. 3(a)], which 
comparable with a typical spectral linewidth of single 




0.0 0.5 1.0 1.5 
Power density (W/cm^; 



DC 
C/D 



b 


Vsi(h) 




r=50K 



- 1 .397 eV 
- 1.412 eV 



-400 -200 200 400 
Rel. exc. energy, E - E^ (|xeV) 



20 40 60 80 100 
Temperature (K) 



FIG. 3: (color online). (a,b) The same as in Fig. 2(a), but 
recorded in the vicinity of Eh with higher spectral resolu- 
tion at different temperatures, (c) The peak-to-peak values 
Ah± of the Vsi(h) ESR lines recorded in magnetic fields Bh- 
and Bh-\- as a function of excitation power density. The solid 
lines are fits to Eq. (2). The dashed lines are linear fits, (d) 
The peak-to-peak values Ah± under excitation into the ZPL 
(Eh = 1.397 eV) and above-ZPL excitation {E = 1.412 eV) 
as a function of temperature. With rising temperature the 
SiC sample becomes impeding the ESR detection in the grey 
shaded temperature region. 



defects [23]. We ascribe these results to different areas 
within our SiC crystal, having slightly varied ZPL ener- 
gies due to different local environment. 

Based on the data of Figs. 2 and 3(a) an important 
conclusion can be drawn: the Vsi(h) spins are only ad- 
dressed when the optical resonance and ESR conditions 
are simultaneously fulfilled. The contrast - i.e., the ra- 
tio of the ESR signal between on and off resonant op- 
tical excitation - is above 200. No signature of Vsi(li) 
ionization is observed. Such a behavior differs from di- 
amonds, where for the resonant control of NV defects 
an additional non-resonant illumination is necessary to 
deshelve the NV defects from the dark state [24] . To the 
best of our knowledge, in the solid state, similar double 
radio-optical resonance has only been observed so far in 
quantum dots [25]. 

We now discuss the efficiency of the optical resonance 
pumping. The Ah± dependence on the excitation power 
density P is presented in Fig. 3c. In the range of P under 
study we observe A^- ~ —Ah-\-^ meaning that the photo- 
induced spin polarization is significantly larger than that 
due to the Boltzmann statistics. The experimental data 
can be well described by a standard model for spin pump- 
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FIG. 4: (color online), (a) ODMR spectra, obtained under 
optical excitation with an energy of the Vsi(h) ZPL tran- 
sition Eh = 1.397 eV and with a slightly detuned energy 
E = Eh + lO/xeV. B\\c. (b) Time-resolved ODMR signal 
(the lower curves) recorded after a RF pulse (the upper curve) 
under the double spin resonance condition, i.e., with an ex- 
citation energy Eh and in a magnetic field Bh-. The thick 
solid lines represent exponential fits Ah exp(— at) for two ex- 
citation power densities P = 1.4 and 0.14 W/cm^. (c) Decay 
rate a as a function of P. A linear extrapolation to P = 
yields the spin-lattice relaxation time Ti = 0.096 d= 0.011 ms. 



ing, which was initially applied for atoms [21, 26] 
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Here, Pq = 2.62 W/cm^ is a characteristic power den- 
sity obtained in time-resolved experiments as described 
later. A phonon- assisted spin pumping with the excita- 
tion energy 15 meV above the ZPL is less efficient and 
for P < IW/cm^ linearly depends on the laser power 
density. 

With rising temperature, the resonance spin pumping 
is observed up to 50 K [Fig. 3(d)]. Remarkably, at T = 
50 K the multiple-line structure is not resolved any more, 
but the contrast still remains the same [Fig. 3(b)]. The 
suppression above 50 K is accompanied by the increase 
of sample conductivity due to the ionization of nitrogen 
donors. The conductivity is probably the only limiting 
factor, as in insulating SiC samples the spin pumping is 
efficient even at room temperature [11]. 

We now discuss possible applications of the double 
radio-optical resonance in SiC. Due to local environment, 
the optical resonance and ESR energies are individual 
for each defect. Additionally, they can be changed by 



a local electric or magnetic field in the range of 50 /ieV 
[23, 27, 28]. This can eventually provide a spectroscopic 
tool to selectively address and manipulate entangled Vsi 
qubits by varying the excitation energy or alternatively 
by tuning the double radio-optical resonance conditions 
for a given Vsi qubit. 

The demonstrative experiments are performed using 
the optically detected magnetic resonance (ODMR) tech- 
nique, which can be sensitive to a single defect spin [29] . 
The ODMR signal is obtained as a normalized change 
in photoluminescence APL/PL of the sideband phonon 
replicas, as indicated in Fig. 1(b). We note that the Vsi- 
related PL is recorded even at room temperature. In 
the ODMR experiments the PL was passed through a 
925 nm longpass filter and detected by a Si photodiode. 
In order to improve the signal-to-noise ratio the RF was 
modulated on/off at a frequency of 4.2 kHz and the pho- 
tovoltage was locked-in. 

Figure 4(a) shows APL/PL as a function of magnetic 
field obtained under the resonance excitation with the 
energy E^. Two Vsi(h) lines are observed at magnetic 
fields Bh- and 5/^+, similar to the light-enhanced ESR 
experiments of Fig. 1(c). When the excitation energy is 
detuned off the optical resonance by only 10 /ieV, the 
ODMR signal is strongly suppressed. 

Finally, we measured spin dynamics of Vsi(h) under 
the double radio-optical resonance conditions, i.e., at E^ 
and Bh- . The experimental details are given in the Sup- 
plemental Material [21]. In brief, the system is optically 
pumped with continuous wave (cw) excitation. An in- 
tense RF pulse equalizes the spin population in the dif- 
ferent sublevels. The time evolution of the ODMR sig- 
nal after the RF pulse follows an exponential decay. Fit 
examples to Ahex.-p{—at) for different pump power den- 
sities P are presented in Fig. 4b. As can be shown [21], 
the decay rate a depends on the spin-lattice relaxation 
time Ti and characteristic power density Pq as 
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A linear extrapolation of a{P) to P = yields Ti = 
0.096 ± 0.011 ms [Fig. 4(c)]. This compares reasonably 
well with that of N V defects in diamonds [30] . 

Conclusions and outlook. — The silicon vacancy defects 
in SiC combine the advantages of semiconductor quan- 
tum dots and the nitrogen-vacancy defects in diamond 
in one system, making them very attractive for quantum 
spintronics applications. They are well isolated from SiC 
lattice, and tunable lasers with narrow linewidth can be 
utilized for the double radio-optical resonance coherent 
control of Vsi qubits, by analogy with atoms. Given that 
silicon and carbon have nuclear spin free isotopes, one 
would expect extremely long spin coherence in isotopi- 
cally purified ^^Si^^C , similar to phosphorous donors 
in pure ^^Si [31]. An important technological aspect of 
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Vsi in SiC is their possibility to be created in state-of- 
the-art transmission electron microscopes (TEMs) [32]. 
Considering that the electron beam in a TEM can be 
focused better than 1 nm, previously unachievable per- 
spectives for spin engineering on an atomic scale may 
become feasible. Furthermore, since spatially separated 
Vsi defects can be selectively addressed even within a 
single SiC nanocrystal, they can be used as a sensitive 
probe to image magnetic field vectors and field gradients 
at nanoscale. And last but not least, single Vsi centers in 
SiC are optically active in the NIR spectral range, over- 
lapping with the first telecom window. This may allow 
for the practical realization of quantum communications 
as silicon vacancies can be incorporated into SiC-based 
LED structures or SiC-based photonic crystals. 
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1 Experimental methods 

1.1 Electron spin resonance (ESR) 
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Figure SI: Sketch of the ESR setup. The investigated sample is put into a microwave resonance 
cavity, surrounded by magnetic coils. Through an optical access a laser beam is coupled into the 
system. The variation of microwave absorption is detected with a lock-in amplifier triggered by a 
small superimposed magnetic field. 
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The experimental setup for performing ESR studies of the SiC samples is sketched in Fig. SI. Via 
a quartz ESR tube, the investigated sample is placed into the center of a microwave resonance cavity 
surrounded by a pair of magnetic coils, which allows for applying a homogeneous external magnetic 
field. The cuboid cavity is the core of the home-modified X-Band spectrometer (Bruker 200D) and 
built in a way, that, if a microwave is coupled into it, standing waves in the transverse electric TE102 
mode will emerge. The indices are related to the number of half-waves in each spatial direction. A 
distinctive feature of this standing wave mode is that its magnetic field maximum oscillates at the 
center of the cavity. The Q-factor (quality factor) of the used cavity is Q = 2800 — 3200. The 
microwave bridge contains the devices for generation and detection of the microwave radiation with 
a frequency of about 9 — 10 GHz, corresponding to A ^ 3 cm. An automated frequency control circuit 
stabilizes the adjusted resonance frequency, which is accurately read out using a frequency counter. 

For the optical excitation of the sample a tunable (875 — 940 nm) diode laser system (Sacher Lion 
TEC 500) is used. The wavelength can be coarsely set with a step motor motion controller, which 
alters the angle of an diffraction grating built into the laser system. The laser cavity can be slightly 
deformed using a piezo element which allows for very fine additional tuning of the laser energy within 
a range of about ±10 GHz. A special feature of the laser beam is the very narrow linewidth of about 
1 MHz, which enables for very accurate investigation of the excitation properties. Via an optical access 
the sample is irradiated with the laser. The beam is widened by a lens mounted in the optical path 
guaranteeing illumination of the whole sample. 
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Figure S2: a, X-Y plot of an ESR spectrum obtained under OR excitation at the Vsi(h) ZPL of 6H-SiC. 
b,c. Corresponding X and Y values of the ESR spectrum from a. Magnetic field is color-coded. 



In order to increase the detection sensitivity the signal obtained by the microwave detector is 
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enhanced by a lock-in amplifier. The external magnetic field generated by the pair of coils is super- 
imposed by a small wobble field with a frequency of 100 kHz. In ESR measurements, the microwave 
absorption is obtained while sweeping the magnetic field. Thus, the modulation results in detection 
of the first derivative of the signal. A certain phase shift 6 of the signal occurs with respect to the 
reference signal of the lock-in, which is ascribed to a time delay caused by the response time of the 
cavity, the electronics and physical procedures. Depending on the phase setting certain amplitudes in 
the ESR spectrum are emphasized, while others are diminished or, in case of their phase being shifted 
by 90°, totally faded out. Hence, a powerful tool to distinguish and identify superimposed signals is 
to analyze the X-Y-plot. This can be seen in Fig. S2a, where the phase diagram of an ESR spectrum 
is presented. The magnetic field in the figure is coded by color, and the signal phases at 332.0 and 
341.0 mT are obviously different from that at 336.0 mT. By variation of the detecting phase it is 
possible to highlight different resonances (Figs. S2b and c). We use this technique to maximize the 
signal amplitude of the outer ESR lines ascribed to the Vsi(h) defect. 

To improve the accuracy of the detection of the external magnetic field a NMR magnetometer is 
mounted next to the cavity. The Larmor frequency of protons in a reference material with known 
g-factor is measured and read out by a frequency counter. 

1.2 Optically detected magnetic resonance (ODMR) 

The ODMR setup shown in Fig. S3 is a variation of the ESR setup presented in the previous section. 
The cavity is designed for optical access from the front and the rear side allowing for coupling of a 
laser beam and optical detection. The photoluminescence (PL) is detected in a transmission geometry 
by a Si photodiode and the laser emission is filtered out by a 925 nm longpass filter. A different 
microwave source is used (Wiltron 69137A Sweep Generator), which features the variation of the 
emitted microwave power and the possibility to chop the microwave by an external driver. The 
maximum accessible output power of the microwave source is 60 mW, which is additionally amplified 
reaching a maximum power within the cavity of 2W. The change in PL intensity induced by the 
change in the ground state population is locked-in. For taking into account that the overall PL is 
governed by different external factors and for enabling to compare different spectra the change in PL 
is normalized, APL/PL . 

In transient ODMR experiments, instead of using a lock-in amplifier the pre-amplified signal is 
detected with an oscilloscope (LeCroy WaveRunner 610zi) triggered by the chopped microwave source. 
In order to improve the signal-to-noise ratio the signal is averaged over 1000 scans. Using a function 
generator the pulse length of the microwave can be varied in order to systematically investigate the 
population difference of the involved spin levels, which is proportional to the obtained ODMR signal 
due to spin-selective PL. 

2 Optical spin pumping of the silicon vacancy 
2.1 Absorption and emission of the silicon vacancy 

Considering a defect center embedded into a crystal, the interaction with light can be explained by the 
interplay of electronic and vibronic transitions between different energy levels (Fig. S4). The transition 
between the lowest vibrational levels of the ground state (GS) and the excited state (ES) not involving 
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Figure S3: Scheme of the ODMR setup. The investigated sample is put into a microwave resonance 
cavity surrounded by magnetic coils. Via optical accesses the sample is illuminated with an infrared 
laser and the emitted light is obtained with a silicon detector. Sweeping the magnetic field the relative 
change in PL intensity is obtained for a chopped microwave irradiation using lock-in detection. The 
transient signals are recorded using an osciloscope. 




Figure S4: Scheme of the resonance excitation (Res. Exc.) and zero-phonon line (ZPL) emission 
between the ground state (GS) and the excited state (ES). Radiative recombination accompanied by 
a photon emission results in phonon replica. The GS has high-spin state with spin multiplicity either 
5=1 or 5 = 3/2. Spin polarization under optical excitation occurs via spin dependent relaxation 
through the metastable state (MS). 
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phonons is called zero-phonon line (ZPL). The phonon transitions of all lattice modes contribute and 
form a continuous phonon replica band of the zero phonon line, which is blue-shifted for absorption 
respectively red-shifted for PL. While a sharp ZPL can be observed at cryogenic temperatures, the 
probability for phonon-assisted transitions increases with temperature due to thermal activation. 

One of the special properties of the investigated silicon vacancy defects constitutes the possibility 
of polarizing their ground state. The working model of the energy structure and optical pumping 
scheme can be adopted from that for NV defects in diamond. The photo-kinetic processes leading 
to spin alignment under optical pumping can be described by the presence of an excited metastable 
state (MS) and a spin-dependent inter-system crossing (ICS) between the high-spin ES, the MS and 
the high-spin GS (Fig. S4). The optical transitions between the GS and the ES are spin-conserving. 
Remarkably, the nonradiative decay via the MS is spin-selective resulting in a polarization of the GS. 
Depending on the polytype of SiC, the lattice site and temperature the transition rates between the 
energy levels vary resulting in different population of the spin sublevels. 

2.2 ESR and ODMR spectra of the silicon vacancy 

The spin Hamiltonian of the silicon vacancy in magnetic field can be written in the form 



Here, Sz is the spin projection on the c axis of SiC crystal. The parameter D characterizes the axial 
crystal field. There is still no agreement in the community whether the ground state of Vsi in SiC has 
spin 5=1 or 5 = 3/2. We now show that the ESR and ODMR spectra presented in this work are 
indistinguishable between these cases. 
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Figure S5: a, A scheme of the light-enhanced ESR in case when the ground state has spin 5=1. b, 
A scheme of the light-enhanced ESR in case when the ground state has spin 5 = 3/2. 
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2.2.1 S = 1 ground state 

Following Eq. (S4), zero field splitting A = D occurs between the rus = and rris = ±1 spin sublevels 
(Fig. S5a). In an external magnetic field B applied parallel to the c axis of SiC crystal the rris = ±1 
sublevel is split further as E±i = ±gefiBB. In the first order, the allowed transitions are with Srris = ±1 
and the ESR for the transitions +1^0 and —1^0 occur at different magnetic fields, labeled in 
Fig. S5a as B- and B^ respectively. There is the following interconnection between these fields and 
A 

hi^ESR = ±A + gel^BBzp . (S5) 

Here, vesr is the resonance frequency of the ESR spectrometer. Under illumination the system 
is pumped into the = ±1 states. This implies inverse population: at B^ the RF radiation is 
absorbed, while at B- the RF radiation is emitted. This is what we observe in our light-enhanced 
ESR experiments (see Figs. Ic and d of the main text). 

Remarkably, the average photoemission rate involving the rris — and rris = ±1 spin sublevels is 
different. Resonant RF changes spin population, resulting in a change of the PL intensity (APL/PL). 
The corresponding ODMR spectrum should reveal two lines, as we observed in our experiments (see 
Fig. 4a). 

2.2.2 S = 3/2 ground state 

Following Eq. (S4), zero field splitting A = 2D occurs between the rris = ±1/2 and rris = ±3/2 spin 
sublevels (Fig. S5b). In an external magnetic field B applied parallel to the c axis of SiC crystal the 
rris = ±1/2 and rris = ±3/2 sublevels are split further as £'±1/2 = ±^e/^5^/2 and £^±3/2 = ±3^eMB^/2 
respectively. In the first order, the allowed transitions are with Srris = ±1 and the ESR for the 
transitions +3/2 +1/2 and —3/2 —1/2 occur at different magnetic fields, labeled in Fig. S5a as 
B- and B^ respectively. There is the following interconnection between these fields and A 

huESR = ±A + QefJ^BB^ , (S5) 

which is equal to Eq. (S5). Under illumination the system is pumped into the rris = ±3/2 states. 
Again, this implies inverse population: at B^ the RF radiation is absorbed, while at B- the RF 
radiation is emitted. The third allowed transition —1/2 +1/2 should appear at the magnetic field 
Bq = [Bj^ + 5_)/2. However, it is insensitive to the illumination and indistinguishable from other 
resonances appearing in ESR spectra at Bq (see Figs. Ic and d). 

The average photoemission rate involving the rris = ±1/2 and rris = ±3/2 spin sublevels is different. 
Therefore, redistribution in the spin population between the rris — —1/2 and rris — +1/2 sublevels 
does not result in a change of the PL intensity and the corresponding ODMR spectrum should again 
reveal two lines (see Fig. 4a). 

2.3 Model of the silicon vacancy resonant spin pumping 

We consider a two- level system with spin populations denoted as n+ and n_. In case of the 3 = 1 
ground state, n+ corresponds to the aggregate population of the rUs = ±1 sublevels and n_ corresponds 
to the population of the = sublevel. In case of the 5 = 3/2 ground state, n+ corresponds to the 
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aggregate population of the = ±3/2 sublevels and n_ corresponds to the aggregate population of 
the rris = ±1/2 sublevels. 

An intense RF pulse equalizes the spin population and at time t = we have n+(0) = n_(0) = ^. 
After the RF pulse the system is optically pumped, i.e., the population n+ increases and the population 
n_ decreases. The spin pumping rate is higher the higher the laser power density P and the higher 
the population n_ are. Simultaneously spin relaxation with the spin-lattice relaxation time Ti occurs. 
It is proportional to the difference from the population in the thermal equilibrium neq±^ which can be 
well approximated to rieq- ^ rieq^ ^ ^ . Finally, the following differential equations can be written 

= kPn. - 1 fn+ - ^ V (S5) 



dt ^ Ti V 2 

"^-n- = -kPn- -^(n--l] . (S5) 



dt Ti V 2 

Here, is a coefficient of proportionality, such that kPri- defines the spin pumping rate. The ESR 
signal is proportional to the population difference An = n+ — n_ , and with n+ (t) + n_ (t) = 1 we 
obtain 

d . , ^ /, ^ 1 



An = kP- [kP ± ^ J An . (S5) 
With boundary condition An(0) = the solution of Eq. (S5) is 

An{t) = Anoo (l - e""^) (S5) 

with spin relaxation rate 

Ti V Po 



a=^(l + ^] (S5) 



and steady-state value 
Here, Pq = ^ is a characteristic power density. 



